Regularizing cosmological singularities by varying physical constants. 
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Varying physical constant cosmologies were claimed to solve standard cosmological problems such 
as the horizon, the flatness and the A-problem. In this paper, we suggest yet another important 
application of these theories: solving the singularity problem. By specifying some examples we show 
that various cosmological singularities may be regularized provided the physical constants evolve in 
time in an appropriate way. 
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I. INTRODUCTION 

One of the most intriguing problems in cosmology is 
the problem of singularities. They are very well-defined 
in relativity and are shown to appear under quite gen- 
eral conditions of geodesic incompletness and a blow-up 
of various geometrical and physical quantities [l| . Up to 
a very recently, the main concern of cosmologists was the 
big-bang singularity in the past which seemed to be un- 
avoidable both in relativity and also in extended scalar 
field theories Q. A lot of generalized theories which in- 
cluded gravity were proposed in order to avoid big-bang. 
Among them the supcrstring and the brane theories [3j,|4|, 
loop quantum gravity ||, higher-order gravity @, and 
many others. The main achievement of such approaches 
was the extension of the evolution of the universe through 
a big-bang singularity like in the pre-big-bang [3] and the 
cyclic § scenarios. 

After the discovery of the accelerated expansion of 
the universe deeper studies of the phenomenon 

of the dark energy showed the plethora of new singu- 
larities ("exotic" singularities) different from big-bang. 
Firstly, a big-rip associated with the phantom dark en- 
ergy was investigated [12j ] , and further classified as type 
I [13| . Then, a sudden future singularity (SFS or type II) 
was proposed [Tij as well as numerous other types such 
as: generalized sudden future singularities (GSFS), finite 
scale factor singularities (FSF or type III) , big-separation 
singularities (BS or type IV) and w-singularities (type V) 
[lot [Tr|. The singularities which fall outside this classifi- 
cation (with perhaps a big-bang as type [16]) are curva- 
ture singularities with respect to a parallelly propagated 
basis (p.p. curvature singularites) which show up as di- 
rectional singularities [17| and also intensively studied 
recently: the little- rip singularities [ijj], and the pseudo- 
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rip singularities [19j . All the above singularities are char- 
acterized by violation of all, some or none of the energy 
conditions which results in a blow-up of all or some of the 
appropriate physical quantities such as: the scale factor, 
the energy density, the pressure, and the barotropic in- 
dex (for a review see Ref. [13] )• In order to be clear, we 
remind that there are three energy conditions: the null 
(gc 2 + p > 0), weak (gc 2 > and gc 2 + p > 0), strong 
(gc 2 + p > and gc 2 + 3p > 0), and dominant energy 
(gc 2 > 0, — gc 2 <p< gc 2 ) (here c is the speed of light, g 
- the mass density in kg m~ 3 , and p - the pressure). 

It emerged fascinating that some of these singularities 
are weaker than big-bang (for example particles [2l| and 
even extended objects [221 ] may pass through them) and 
may appear in the very near future [23l . |24j . Then, it 
is interesting to discuss if there are any physical reasons 
which can "weaken" or just remove a big-bang (or other - 
"stronger" ) singularity from the evolution of the universe. 
Our suggestion here is to make use of the time-evolution 
of the physical constants combined with the dynamical 
evolution of particular models of the universe. It has 
been recently shown [25[ that the quantum effects may 
change the nature of the sudden future singularity as well 
as the big-rip and the finite scale factor singularity. It 
is then reasonable to think of the time-varying physical 
constants to do the job as well. 

The idea of variation of physical constants has been 
established widely in physics both theoretically and ex- 
perimentally. From the theoretical side the early ideas 
of Weyl [26| and Eddington 27| were most successfully 
followed by Dirac's Large Number Hypothesis 28] from 
which it was concluded that the gravitational constant 
should change in time as G(t) cx l/t. This led to 
the scalar-tensor gravity theory developed by Brans and 
Dicke |29| who followed the ideas of Jordan [30|. These 
ideas were further embedded into superstring theories in 
which the coupling constant of gravity became running 
during the evolution of the early universe [3]. In fact, a 
lot of physical constants such as the gravitational con- 
stant G, the charge of the electron e, the velocity of 
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light c, the proton to electron mass ratio rn n / rn P . and 
the fine structure constant a are interrelated |3ll . |32| and 
the variation of one of them may be associated with vari- 
ation of others. However, apart from Brans-Dicke type of 
gravitational constant evolution models, the most popu- 
lar theories which admit physical constants variation are 
the varying speed of light theories (33j and varying alpha 
(fine structure a) theories [34j . It has been shown that 
both of these theories allow the solution of the standard 
cosmological problems such as the horizon problem, flat- 
ness problem, and the A— problem. Here, we will apply 
these theories to solve yet another problem - the singu- 
larity problem. 

The paper is organized as follows. In Section |TT] we 
discuss the basics of varying constants models using the 
new form of the scale factor which encompasses quite 
a few types of singularities after a specific choice of its 
parametrization. In Section IIIII we show the examples of 
regularization the singularities due to the time-variation 
of physical constants. In Section llVI we give our conclu- 
sions. 



II. VARYING CONSTANTS MODELS 

Following the Refs. [HI, we consider the Friedmann 
universes within the framework of varying speed of light 
theories (VSL) and varying gravitational constant theo- 
ries. The field equations read as 



g(t) = 
P(t) = 



kc 2 (t) 



8nG(t) \a 2 
c 2 (t) 



a a 2 kc 2 (t) 
8nG(t) \a + ~^ + a 2 



(11.1) 
(11.2) 



and the energy-momentum conservation law is 



G(t) kc(t)c(t) 
-Q(t)-p7TX + 3- 



G(t) 



A-kGo 2 



(II.3) 

Here a = a(t) is the scale factor, the dot means the 
derivative with respect to time t, G — G(t) is time- 
varying gravitational constant, c = c(t) is time- varying 
speed of light, and the curvature index k — 0, ±1. In 
Ref. [33| the barotropic equation of state for matter was 
assumed. Since we want to discuss more general cases in 
order to obtain various types of singularities, we will not 
assume any link between the energy density and pressure 
(i.e. the equation of state). 

In contrast to many references dealing with sudden fu- 
ture singularities [13, HH , which consider the scale factor 



a(t) 
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, (II.4) 



with S,t s ,a s ,m,n being constants, we propose a new 
form of the scale factor, which after appropriate choice 
of parameters admits big-bang, big-rip, sudden future, 



finite scale factor and ^-singularities and reads as 
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(II.5) 



with the constants t s , a s , m, n. Its expansion around t \ 
t s gives (cf. HI) 



a(t) = a s — nia s 
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Notice that in equation (111.51) , a s has the unit of length 
and all the terms are dimensionless. The scale factor is 
zero (a = 0) at t = (a big-bang singularity), and it 
is a constant a — a s at t — t s . The first and second 
derivatives of the scale factor 5[) are 



a(t) = a(t) 
&(t) = a(t) 
+ a(t) 
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From (|IL7|) - (|II.8|> . one can see that for 1 < n < 2 d(0) = 
oo and a(t s ) — ma s /t s =const., while a(t s ) = a s , a(0) = 
co and d(t s ) = — oo (p — > oo) and we have a sudden 
future sing ularity [Hj]. 

For a flat (k = 0) Friedmann model it is possible 
to write down an explicit relation between the pressure 
and the energy density, though with a time-dependent 
barotropic index, in the form 



p s (t) = w s (t)g s (t) 



where 



w s (t) 



Kt) 



[2q(t) - 1] 



(II.9) 



(11.10) 



and q(t) = —aa/a 2 is the deceleration parameter. 

Using pTT|) . dlL2l . (flTT7|) . and pL8)l for the curvature 
index k = we have 



g(t) 
Pit) 



8irG(t) 
c 2 (t) 
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t ts \ ts 
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(11.11) 
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ts 



n-2 



From (|II.5|) . pi. lit and (III.12|) we can conclude that 
for < m < 2/3 we deal with a big-bang singularity and 
a— >0, g — > oo, p — >ooatt— >0, while for m < we 
have a big-rip singularity with a— >-oo,£>— >-oo,p— >-ooat 
t = 0. Besides, it is clear that a sudden future singularity 
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FIG. 1. The plots of the scale factor a(t), the energy density g(t), and the pressure p(t) for the two specific models given by the 
scale factor (|II.5|) . The first model is for the parameters m — 0.6, n = 1.5 and describes the sudden future singularity (SFS) 
scenario. The second model is for the parameters m — 0.6 and n = 0.5 and describes the finite scale factor singularity (FSF) 
scenario. 



(SFS) which appears for 1 < n < 2 at t = t s (a = a s , g = 
const., p — > oo) and a stronger [T3, HH [33 finite scale 
factor singularity (FSF) appears for < n < 1 at t = t s 
(a = a s , g — > 00, p — > 00). The example plots of these 
models are given in Fig. [1] In fact, for the former only the 
last term in the pressure of the type (1 — t/t s ) n ~ 2 blows- 
up, while for the latter two more terms (1 — t/ts) 71-1 and 
(l-t/t s ) 2 ( n ~V in (|TLT2|) do [35]. Su ch a c hoice also leads 
to a blow-up of the energy density pi. Ill) . 

For the sake of further discussion, we will split the scale 
factor (|II.5|) into the two factors, one giving a big-bang 
singularity as b , and another giving an exotic singularity 
a ex , as follows: 



a(t) = a B B ■ a ex , 



where 



a BB 



a ex = a s exp ( 1 - — 



(11.13) 

(11.14) 
(11.15) 



From (III.13|) one sees that in a special case n = 0, one 
recovers a standard big-bang scale factor with 



QBB{t) 

Pbb (t) 



3 m 2 
8ttG(<) 72" ' 

c 2 (t) m(3m-2) 
~8irG(t) t 2 



(11.16) 
(11.17) 



and it is possible to write down an equation of state in 
the form of a barotropic perfect fluid as 

Pbb = (7 - l)Qc 2 (t) = o fJtL , (H.18) 



where 



7 



8TTG(t)t 2 37 



3m 



(11.19) 



as in the standard notation [38|. The standard big-bang 
models given by 14|) are decelerating for 7 > 2/3 (m < 
1), and accelerating for 7 < 2/3 (m > 1). The pressure is 
positive for 7 > 1 and negative for 7 < 1. This of course 
also refers to its value at a big-bang singularity (t = 0), 
where it is plus infinity for 7 > 1, and minus infinity for 
7 < 1. 

Now, let us notice that in the limit m = the Eq. (|II.5P 
reduces to an exotic singularity scale factor given by 
pTT5|) . From pLTTT) and pLF2l) we have 



g e x(t) 

Pex(t) 



87rG(t) t 2 
c 2 (t) 
8nG(t) ' 

.,,2 



2(n-l) 



(11.20) 



(11.21) 
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and so 



/,\ Pex(t) 2/,\ 
Wexit) = — - = -C It) 



Qex{t) 



2n-l 



3 n 
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with the deceleration parameter equals to 



Qex(t) = -1 - 



1 



1 



(11.22) 



(11.23) 



From pi.20jl - pi.22p we may conclude that for rc > 2 the 
energy density and pressure vanish (g = 0, p = 0), while 
the w-index blows-up to infinit y w hich is exactly the char- 
acteristics of a w-singularity 15|. In fact, the tu-index 
blows-up for any positive value of n > (n ^ 1), which 
together with the fact that the energy density diverges 
for < n < 1, and both the energy density and pressure 
diverge for 1 < n < 2, means that we have neither FSF 
nor SFS singularity here. 

Another way of writing clbb is the way one does in 
superstring cosmology Q 



a(t) 



t 



(11.24) 



in which there are four branches (two of them "pre-big- 
bang"). 



which is a faster decrease that in the standard Dirac's 
case G(t) oc \jt [H, HH, would presumably remove a 
big-bang singularity in Friedmann models (i.e. removes 
both p and g singularities). Such a time-dependence of 
G would perhaps be less influenced by the geophysical 
constraints on the temperature of the Earth as it was 
discussed early in Ref. 39]. On the other hand, in the 
Dirac's case G(t) cx 1/t, only the g singularity can be 
removed. 

Another suggestion is that the scale factor PL5J) would 
not approach zero at t — > if it was rescaled be a "regu- 
larizing" factor a rg = (1 + \/t rn ) (m > 0), i.e., 



1 



1 



1 



(IIL2) 



So far we have preliminary discussed a possibility for 
the gravitational constant to evolve in time. Now, we will 
also discuss the time-evolution of the speed of light. It is 
clear from pi.lip - pi.12l) that any change of the type of a 
singularity which is also singularity of density cannot be 
done without admitting a curvature term in the Einstein 
equations (|II.1|) - (|II.2|> . This, especially refers to a big- 
bang singularity - it cannot be removed at all, unless the 
spatial curvature is non-zero. 

b) other exotic singularities 

In order to regularize an SFS singularity by varying 
speed of light we suggest that the time-dependence of 
the speed of light is given by 



III. REGULARIZING SINGULARITIES 

Our main idea is to investigate how different types of 
time-evolution of the physical constants like c and G in- 
fluence the evolution of the universe. In particular, we 
are interested in a possibility to change the nature of sin- 
gularities due to the variability of these constants. The 
review of the experimental bounds on the variability of 



the constants are given in many references - the Ref. [31 1 
is one of the most recent ones. 

It is quite reasonable that the nature of singularities 
may change while the constants are evolving. In fact, an 
analogous phenomenon, though due to a different physi- 
cal reason - the conformal anomaly - results in strength- 
ening a singularity |25|. Namely, an SFS singularity be- 
comes an FSF sin gula rity - the latter is a stronger type 
of singularity [l7l l2ll. l3a |. However, in our further in- 
vestigations we will be discussing mostly weakening the 
singularities - this is a natural expectation in cosmology 
which helps to solve the singularity problem. 

a) a big-bang singularity 

For the purpose described above, we first suggest that 
the time- variation of the gravitational constant in pi. lip 
and pTT2|) of the form 



G{t) oc 



t 2 



(111.1) 



c(t) 



(III.3) 



(cq = const., p = const.) which after substituting into 
(IILT2P gives 
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(III.4) 



It then follows from (|III.4[) that an SFS singularity is 
regularized by varying speed of light provided that 



p>2-n (1 < n < 2) 



(III.5) 



However, there is an interesting physical consequence of 
the functional dependence of the speed of light pil.3l) . 
Namely, it gradually diminishes reaching zero at the sin- 
gularity. In other words, the light slows and eventually 
stops at an SFS singularity. Such an effect was discovered 
in a dense fluid in solid state physics some years ago. 

Note that the regularization of a w-singularity (|II.22|1 - 
PI.23P is possible only by assuming a time- varying speed 
of light pTL3]) only, if 



p > n 



(III.6) 
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One of the standard assumptions on the variation of 
the speed o f lig ht is that it follows the evolution of the 
scale factor [3J 



and an FSF singularity (0 < 1 < n) is regularized when 
r > l-n . (111.16) 



c(t) = c Q a s {t) 



(III.7) 



with Co and s constant. The field equations and 
(111.21) become 



g(t) 



8irG(t) \a 2 

2s 



+ kcW {s - 1} 



(III.. 



p(*) = - C » 



87rG(t) 



2- + ^ + fccoV^ 1 )) .(III.9) 



With such a dependence it is possible to remove a pres- 
sure singularity provided s > 1/m — 3j/2 for k = and 
m > 0, s > 1/2 or m < 0, s < 1/2 for fc 7^ 0. A more 
sophisticated choice of the time variation of the speed of 
light which could regularize the singularities is 



c(t) = co 1 



t 



a s (t) 



or in terms of the scale factor 



c 



i-i 

t. 



(111.10) 



(III.ll) 



Since g(t) does not depend on c{t) (for k = 0), then it 
is impossible to strengthen an SFS singularity to become 
an FSF singularity. It is possible only, if we assume that 
the gravitational constant G changes in time. Let us then 
assume that 



(7(/, = (y„(i-i 



(111.12) 



(r = const., Go = const.) which changes (|II.11[) and 
([Tr[2"]) to the form 



g(t) = 



8ttG 
2 



777 



1 - 



r+2n-2 



2mn 

~tt7 



1 - 



p(t) 



6- 



87rGo 

777,77 



777(3777 — 2) 

r+? 1 — 1 



+ 2 



77(77 — 1) 



i 

1 



1 

77 2 , 
- 3 ^' [ 



r+n-2 



r+n— 1 



(111.13) 



r+2n-2 



(III. 14) 



From dnrm and (IIII.14I) it follows that an SFS singu- 
larity (1 < 77 < 2) is regularized by varying gravitational 
constant when 



r > 2 — 77 



(111.15) 



On the other hand, assuming that we have an SFS sin- 
gularity and that 



1 < r < , 



(111.17) 



we get that varying G may change an SFS singularity 
onto a stronger FSF singularity when 



< r + 77 < 1 



(111.18) 



A physical consequence of the functional dependence of 
the gravitational constant in 1 2[) is that the strength 
of gravity becomes infinite at the singularity. This is 
quite reasonable if we want to regularize an infinite 
(anti)tidal force at the singularity |3a |. 

A hybrid case which would influence both types of sin- 
gularities (big-bang (cf. and other exotic ones) 
is 



G(t) 



Go 
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t 



t 2 \' t s 

which changes ([II. Ill) and (|II.12[) into 
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(111.19) 



r+n— 1 



(111.20) 



r+2n-2 



(111.21) 



From pil.201) and (|III.21[) one can see that a big-bang 
singularity at t = is regularized with no additional 
conditions while SFS or FSF singularities are regularized 
under the conditions (|III.15|) and (|III.16[) appropriately. 

c) singularities in (anti-)Chaplygin gas 
cosmology 

Since a couple of years, there has been a proposal that 
the dark energy can be simulated by a Chaplygin or an 
anti-Chaplygin gas model. One of the interests is that an 
anti-Chaplygin gas model allows the so-called big-brake 
singularity [36] (g — > and p — > 00) which is a special 
case of a sudden future singularity g —¥ const, and p — > 
00. The equation of state of the (anti-) Chaplygin gas 
reads as 



(111.22) 
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where the "-" sign refers to a Chaplygin gas while the "+" 
sign refers to an anti- Chaplygin gas case. It is interesting 
to note that the only limit of an SFS to get a big-brake 
for the scale factor 5[) bearing in mind the fact that 
g — > (i.e., a — > - see pi.7[l ) would be m — > and n — > 
which is not physically interesting. The reasonable limit 
of the scale factor (|II.4[) is possible [23,'}. 
Inserting (IIII.22I) into p!3)) gives [33] 



g(t) 



a V g{t)c 2 (t) 



G(t) , q kc(t)c(t) 
G(t) 



4irG(t)a 2 ' 
(111.23) 

In order to find an exact solution, we will first assume 
that the speed of light is constant (c = const.), and that 
the gravitational constant changes as 



G(t) = G g(t) , 

which gives (|III.23|) in the form 

2c 2 gg = _ 3 a 
g 2 c 2 =F A a 



(111.24) 



(111.25) 



and integrates to give (we should make an assumption 
A < g 2 c 2 for "-" sign which corresponds to a Chaplygin 
gas) 



g(t) 



±A C + 



lTA c 
aHt) 



(111.26) 



with go — const, and A c = A/ go. From (|III.26|) one 
can see that the standard big-bang limit g ~ a~ 3 / 2 is 
achieved if A c — > and A c — > ±1 gives the cosmological 
constant limit. From (IIII.22[) and pil.261) we have 



p(t) = Tc- 



goA c 



±A n 



(111.27) 



±1 



The standard big-bang limit p is given for A c 
and p = ^fc^/go for A c ±1. 

The more interesting solution in order to demonstrate 
regularization of singularities can be obtained in a more 
general case of both varying G — G(t) and c = c(t) 
though with zero curvature k = case by making the 
following ansatz in (jIII.23[) 



g(t)c(t) = B = const. 



which gives 



Q +3 a(B 2 T A 



^ = 
g a \ B 2 ) G 

The solution of (|III.29|) reads as 

g(t)a 3 ~> (t)G(t) = E = const 

where we have defined 

B 2 tA 
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(111.28) 



(111.29) 



(111.30) 



(111.31) 



for the sake of comparison with the standard cosmol- 
ogy (here 7 has not the meaning of the barotropic index, 
but as we can learn it makes the scale factor scaling the 
same as in standard case). Bearing in mind pil.28[) and 
(|III.30[) we obtain that 



c(t) = ~a*i(t)G® . 
From pll.22|) and pll.32|) we finally have 
p(t)=T^c(t) = T^(t)G(t) 



(111.32) 



B 



E 



c(t) a 3 i{t)G(t) ' 
Putting the standard big-bang scale factor pl,14| 



(111.33) 
(111.34) 



a(t) 



t \ 3 ^ 



J \"S 

we now have from (|III.33|) and (|III.34|) 

Et 2 



P(t) 



t 2 G(t) ' 



G(t) 



(111.35) 

(111.36) 
(111.37) 



which give g —¥ 00 and p(0) = provided G(0) = const. 
7^ 0. The singularity at t = in g can be regularized by 
taking G(t) cx 1/t 2 as in pTLT|) at the expense of having 
a constant pressure (cosmological term) instead of zero 
pressure. 



IV. CONCLUSIONS 

We have shown by specifying some examples that it is 
possible to regularize (avoid) cosmological singularities 
due to variation of the physical constants. Although it 
seems to work in a general framework of varying con- 
stants theories, we have considered this phenomenon in 
the theories with varying speed of light (VSL) c(t), and 
with varying gravitational constant G = G(t). 

For example, in order to regularize a big-bang singular- 
ity, the simple modification of Dirac's relation (G cx 1/t 2 
instead of G oc 1/t) is required. In order to regularize a 
sudden future singularity (SFS), finite scale factor singu- 
larity (FSF) or a w— singularity, some more complicated 
time-dependence for c(t) and G(t) is necessary. We have 
found such a dependence and it was appropriate to a 
newly introduced scale factor given by (III. 51) . Interest- 
ingly, in order to regularize an SFS by varying c(t), the 
light should stop propagating at a singularity. On the 
other hand, to regularize an SFS by varying gravitational 
constant - the strength of gravity has to become infinite 
at a singularity which is quite reasonable because of the 
requirement to overcome an infinite (anti-)tidal forces at 
the singularity. 
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We have also studied the regularization of singularities 
by varying c(t) and G(t) within the framework of (anti-) 
Chaplygin gas cosmology and have shown thet regular- 
ization is also possible in these theories. 

Finally, we hope that variation of the physical con- 
stants which leads to regularization of singularities may 
be useful in the discussions of the multiverse concept giv- 
ing the link through a kind of "fake" singularities to vari- 
ous parts of the universe with different physics. Of course 
our discussion is preliminary and should be continued by 



using appropriate mathematical formalism of both gen- 
eral relativity and particle physics. 

V. ACKNOWLEDGEMENTS 

The authors acknowledge the discussions with Adam 
Balcerzak and Tomasz Dcnkiewicz. This work was sup- 
ported by the National Science Center grant No N N202 
3269 40. 



[1] S.W. Hawking and G.F.R. Ellis, The large-scale struc- 
ture of space-time (Cambridge Univ. Press), Cambridge 
(1999) . 

[2] A. Borde and A. Vilenkin, Phys. Rev. Lett. 72, 3305 
(1994). 

[3] J. Polchinski, Rev. Mod. Phys. 68 (1996), 1245; J. 
Polchinski, String Theory (Cambridge University Press, 
Cambridge), 1998. 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett., 83 (1999), 
3370; Phys. Rev. Lett., 83 (1999), 4690; N. Arkani- 
Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. D59, 
086004 (1999); T. Shiromizu, K. Maeda, and M. Sasaki, 
Phys. Rev. D62, 024012 (2000); M. Sasaki, T. Shiromizu 
and K. Maeda, Phys. Rev. D62, 024008 (2000). 

[5] M. Bojowald, Living Reviews in Rel. 8, 11 (2005). 

[6] T. Sotiriou and V. Faraoni, Rev. Mod. Phys. 82, 451 
(2010); S. Nojiri and S.D. Odintsov, Phys. Rept. 505, 
59 (2011); T. Biswas, T. Koivisto, and A. Mazumdar, J. 
Cosmol. Astropart. Phys. 008, 1011 (2010). 

[7] J.E. Lidsey, D.W. Wands, and E.J. Copeland, Phys. 
Rept. 337, 343 (2000); Gasperini M. and Veneziano G., 
Phys. Rept. 373, 1 (2003). 

[8] Khoury J., Steinhardt P.J., Turok N., Phys. Rev. Lett. 
91 (2003), 161301; ibidem 92 (2004), 031302; Khoury J., 
Ovrut B.A., Steinhardt P.J., Turok N., Phys. Rev. D64 
(2001), 123522; ibidem D66 (2002), 046005. 

[9] S. Perlmutter et al., Astroph. J. 517, (1999) 565; A. G. 
Riess et al., Astron. J. 116, 1009 (1998); A.G. Riess et 
al., Astroph. J. 560, 49 (2001). 
[10] J.L. Tonry et al., Astroph. J. 594, 1 (2003); M. Tegmark 
et al., Phys. Rev. D69, 103501 (2004); R.A. Knop et al., 
Astrophys. J. 598, 102 (2003). 
[11] M. Kowalski et al., Astrophys. J. 686, 749 (2008); Aman- 
ullah et al. (The Supernova Cosmology Project), As- 
troph. J. 716, 712 (2010). 
[12] R.R. Caldwell, Phys. Lett. B 545, 23 (2002); M.P 
Dabrowski, T. Stachowiak and M. Szydlowski, Phys. 
Rev. D 68, 103519 (2003); R.R. Caldwell, M. 
Kamionkowski, and N.N. Weinberg, Phys. Rev. Lett. 91, 
071301 (2003); PH. Frampton, Phys. Lett. B 562 (2003), 
139; H. Stefancic, Phys. Lett. B586, 5 (2004); S. Nojiri 
and S.D. Odintsov, Phys. Lett. B595, 1 (2004). 
[13] S. Nojiri, S.D. Odintsov and S. Tsujikawa, Phys. Rev. D 
71,063004 (2005); K. Bamba, S. Capozziello, S. Nojiri, 
and S.D. Odintsov, arXiv: 1205.3421. 
[14] J.D. Barrow, G.J. Galloway and F.J. Tipler, Mon. Not. 
Roy. Astr. Soc, 223, 835 (1986); J.D. Barrow, Class. 
Quantum Grav. 21, L79 (2004); J.D. Barrow and Ch. 



Tsagas, Class. Quantum Grav. 22, 1563 (2005); V. Sahni 

and Yu.V. Shtanov, Class.Quant.Grav. 19, L101 (2002). 
[15] M. P. Dabrowski and T. Denkiewicz, Phys. Rev. D 79, 

063521 (2009); L. Fernandez- Jambrina, Phys. Rev. D82, 

124004 (2010). 
[16] C. Kiefer, Ann. Phys. (Berlin) 19, 211 (2010). 
[17] L. Fernandez- Jambrina, Phys. Lett. B656, 9 (2007). 
[18] P.H. Frampton, K.J. Ludwick, R.J. Scherrer, Phys. Rev. 

D84, 063003 (2011). 
[19] P.H. Frampton, K.J. Ludwick, R.J. Scherrer, arXiv: 

1112.2964. 

[20] M.P. Dabrowski and T. Denkiewicz, AIP Conference Pro- 
ceedings 1241, 561 (2010); arXiv: 0910.0023. 

[21] L. Fernandez- Jambrina and R. Lazkoz, Phys. Rev. 
D70, 121503(R) (2004); L. Fernandez- Jambrina and R. 
Lazkoz, Phys. Rev. D74, 064030 (2006). 

[22] M.P. Dabrowski and A. Balcerzak, Phys. Rev. D73, 
101301(R) (2006). 

[23] M.P. Dabrowski, T. Denkiewicz, M.A. Hendry, Phys. 
Rev. D75, 123524 (2007); H. Ghodsi, M. A. Hendry, M. 
P. Dabrowski, T. Denkiewicz, Mon. Not. Roy. Astron. 
Soc. 414, 1517 (2011); T. Denkiewicz, M.P. Dabrowski, 
H.G. Ghodsi, and M.A. Hendry, Phys. Rev. D85, 083527 
(2012). 

[24] T. Den kiewicz, J. Cosm ol. Astropart. Phys. (2012), to 
appear; larXiv:1112.5447l 

[25] M.J.S. Houndjo, Europhys. Lett. 92, 10004 (2010); 
M.J.S. Houndjo et al. arXiv: 1203.6084. 

[26] H. Weyl, Ann. Phys., 59, 129 (1919); Natturwis- 
senschaften 22, 145 (1934). 

[27] A.S. Eddington, The Mathematical Theory of Relativity, 
(Cambridge University Press, Cambridge), 1923; New 
Pathways in Science, (Cambridge University Press, Cam- 
bridge), 1934. 

[28] P. A.M. Dirac, Nature 139, 323 (1937); Proc. Roy. Soc. 

A165, 189 (1938). 
[29] Brans, C, and Dicke, R.H., Phys. Rev. 124 (1961), 925 . 
[30] P. Jordan, Zeit. Phys. 157, 112 (1959) . 
[31] J.-P. Uzan, Rev. Mod. Phys. 75, 403 (2003); arXiv: 

1009.5514 . 

[32] J.D. Barrow, The Constants of Nature, (Vintage Books, 
London), 2002. 

[33] J.D. Barrow, Phys. Rev. D59, 043515 (1999); A. Albrecht 
and J. Magueijo, Phys. Rev. D59, 043516 (1999); J.D. 
Barrow and J. Magueijo, Class. Quantum Grav. 16, 1435 
(1999); P. Gopakumar and G.V. Vijayagovindan, Mod. 
Phys. Lett. A16, 957 (2001). 



8 



[34] J.D. Barrow and J. Magueijo, Phys. Lett. B443, 104 
(1998); J.D. Barrow, Ann. Phys. (Berlin), 19, 202 (2010). 

[35] F. Tipler, Phys. Lett. A 64, 8 (1977); A. Krolak, Class. 
Quantum Grav. 3, 267 (1988). 

[36] V. Gorini, A. Yu. Kamcnshchik, U. Moschella, and V. 
Pasquier, Phys. Rev. D69, 123512 (2004); Z. Li, P. Wu, 



and H. Yu, J. Cosmol. Astropart. Phys. 09, 017 (2009). 
[37] A.K. Singha and U. Debnath, Int. J. Mod. Phys. D16, 
117 (2007). 

[38] M.P Dabrowski, Phys. Lett. B702, 320 (2011). 
[39] E. Teller, Phys. Rev. D73, 801 (1948). 



